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ABSTRACT 

We have carried out the first systematic survey of the sub-millimeter properties of broad absorption 
line (BAL) quasars. 30 BAL quasars drawn from a homogeneously selected sample from the Sloan 
Digital Sky Survey at redshifts 2 < z < 2.6 were observed with the SCUBA array at the JCMT 
to a typical rms sensitivity of 2.5 mJy. Eight quasars were detected at > 2a significance, four of 
which are at > 3cr significance. The far-infrared luminosities of these quasars are > Lq. There 
is no correlation of sub-millimeter flux with either the strength of the broad absorption feature or 
with absolute magnitude in our sample. We compare the sub-millimeter flux distribution of the BAL 
quasar sample with that of a sample of quasars which do not show BAL features in their optical 
spectra and find that the two are indistinguishable. BAL quasars do not have higher sub-millimeter 
luminosities than non-BAL quasars. These findings are consistent with the hypothesis that all quasars 
would contain a BAL if viewed along a certain line-of-sight. The data are inconsistent with a model 
in which the BAL phenomenon indicates a special evolutionary stage which co-incides with a large 
dust mass in the host galaxy and a high sub-millimeter luminosity. Our work provides constraints on 
alternative evolutionary explanations of BAL quasars. 

Subject headings: galaxies: active - galaxies: evolution - quasars: absorption lines - submillimeter 



1. INTRODUCTION 

The striking correlation between black hole mass and 
the mass of stars in spheroidal galaxies suggests a physi- 
cal link between the build up of black holes by accretion 
and the formation of stars. Most models for this link 
(e.g. Silk & Rees 1998; Fabian 1999) invoke feedback 
mechanisms in which black holes drive outflows: once 
the black hole reaches a critical mass fixed by the host 
galaxy mass, these outflows drive gas and dust out of 
the galaxy and terminate both star formation and fur- 
ther black hole growth. 

Observational evidence for outflows in quasars comes 
in the form of broad absorption lines (BALs). These are 
metal absorption systems with huge blue-shifted veloci- 
ties (several thousand kms^^). The similarity of the con- 
tinuum and line emission of BAL and non-BAL quasars 
motivates the hypothesis that BAL quasars are not in- 
trinsically different from other quasars. The presence 
of BAL features in the spectra of only 15% of optically 
selected quasars (Tolea et al. 2002) could naturally be 
explained by a difference in viewing angle if the sub- 
relativistic outflow at the origin of the BAL features is 
not isotropic (Weymann et al. 1991). The popular no- 
tion that the outflow is preferentially located close to 
the edge of the torus surrounding the supermassive black 
hole has found support in spectropolarimetric measure- 
ments (Goodrich & Miller 1995; Cohen et al. 1995). 
However, studies of radio-loud BAL quasars, in which 
the radio properties give some information on the orien- 
tation, appear at odds with this simple orientation model 
(Becker et al. 2000). The other main contender as an 



interpretation of the BAL phenomenon is the evolution- 
ary scenario of Briggs et al. (1984) in which all quasars 
pass through a BAL phase for ~ 15% of their active life- 
times. The small sizes of the radio lobes in radio-loud 
BAL quasars (Becker et al. 2000) is suggestive of the 
BAL phase co-inciding with an early stage of quasar ac- 
tivity - perhaps removing a shroud of gas and dust from 
the nuclear region. 

The sub-millimeter emission from quasars comes from 
optically-thin, cool dust and is therefore orientation- 
independent. It is also likely heated by young stars 
in starbursts which will evolve over the lifetime of the 
quasar. It is therefore the prime wavelength in which 
to discriminate between the orientation and evolutionary 
explanations of the BAL phenomenon. If all quasars con- 
tain BALs, then BAL quasars would be expected to have 
the same sub-millimeter properties as non-BAL quasars. 
But if the BAL phenomenon is a phase that all quasars 
go through, and is connected with the termination of 
large-scale star-formation, then the sub-millimeter emis- 
sion should differ between BAL and non-BAL quasars. 

There are hints from previous studies with millime- 
ter and sub-millimeter photometry that BAL quasars 
are more luminous than non-BAL quasars, but the num- 
bers in these studies are small. Omont et al. (1996) 
found that 2 out of 4 BAL quasars were detected at 1.1 
mm, whereas only 4 out of 19 non-BAL quasars were 
detected. More recently, Carilli et al. (2001) observed 
41 quasars from the Sloan Digital Sky Survey (SDSS) 
with IRAM. However, their sample contained only 7 BAL 
quasars with a detection fraction of 4/7. Compared with 
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their non-BAL detection fraction of 12/34, this is con- 
sistent with BALs being more s\ib-millimeter luminous, 
but the numbers are too small to provide a significant 
result. Priddey et al. (2003) detected 2 out of 8 BAL 
quasars with SCUBA at the JCMT, compared to a de- 
tection rate of 7 out of 49 non-BALs - again suggestive 
of BALs being more sub-millimeter luminous, but not 
at a significant level. Willott et al. (2002) and Rawl- 
ings et al. (2003) both found that the only moderate 
radio-power quasars in their samples with luminous sub- 
millimeter emission were the sole objects in each sample 
to possess BALs. 

The Sloan Digital Sky Survey (SDSS) is transform- 
ing the study of quasars with the size, data quality and 
completeness of its quasar survey. It is now possible to 
define very large samples of BAL quasars with homoge- 
neoiis selection criteria. BAL quasar samples have been 
derived from the SDSS Early Data Release (Stoughton 
et al. 2002) by Tolea et al. (2002) and Rcichard et al. 
(2003). In this paper we present sub-millimeter obser- 
vations of a carefully selected sample of broad absorp- 
tion line quasars from the SDSS. These arc compared 
to quasars which do not possess BALS with similar red- 
shift and luminosity to determine if there is a difference 
between the sub-millimeter properties of BAL and non- 
BAL quasars. 

We assume throughout that Ha = 70kms~^ Mpc^^, 
Hm = 0.3 and ^Ia ~ 0.7. All data used from papers in 
the literature have been adjusted to be consistent with 
this cosmological model. In Sec. 2 we describe the sub- 
millimeter observations performed on the sample of BAL 
quasars. In Sec. 3 we compare the sub-millimeter prop- 
erties of BAL quasars with those from a matched sample 
of quasars without BALs. In Sec. 4 we interpret these 
results and present our conclusions in Sec. 5. 

2. SUB-MILLIMETER PHOTOMETRY OF BAL QUASARS 

2.1. BAL quasar sample selection 

In order to determine if BAL quasars have different 
sub-millimeter properties from non-BAL quasars, we re- 
quire a control sample of non-BAL quasars with matched 
luminosity and redshift. To avoid doubling our required 
observing time by defining our own control sample, we 
instead defined our BAL sample so that it is similar to 
non-BAL sub- millimeter /millimeter quasar surveys per- 
formed by others. Since BAL quasars are usually defined 
on the basis of absorption of the C IV or Mg li lines, highly 
complete samples are limited to z < 4 (e.g. Reichard et 
al. 2003). Further, since the SDSS contains many more 
quasars at 2; ~ 2 than ^ ~ 3, to define a large sample 
of BAL quasars we need to work at redshifts ~ 2. Two 
recent quasar surveys have been carried out at this red- 
shift with SCUBA at the JCMT (Priddey et al. 2003) 
and MAMBO at the IRAM 30-m (Omont et al. 2003) 
and we combine these for the control sample (further de- 
tails are given in Sec. 3.1). 

To select a sample of BAL quasars for observation we 
use the SDSS EDR BAL sample of Reichard et al. (2003) , 
in preference to that of Tolea et al. (2002), since it uti- 
lized a more detailed selection method. Our BAL quasar 
sample is restricted to the redshift range 2 < z < 2.6 
to be comparable with the bulk of sources in the control 
sample. We further restrict the sample to the RA range 
10^ - 18^ for compatibility with telescope scheduling. 



The balnicity index Bl gives a measure of the strength of 
broad absorption features (e.g. Weymann et al. 1991). 
All the quasars in our redshift range show the C IV wave- 
length region in their spectra, so our selection is on the 
basis of the Civ balnicity. The Reichard et al. sample 
contains 39 quasars in this redshift and RA range with a 
wide range of BI (0 — 10 000 km s^^). The four quasars 
with BI < 200kms~^ were excluded from our sample 
since these represent very weak broad absorption fea- 
tures. Of the remaining 35 quasars we selected a sampk; 
of 30 with the highest dercddened absolute magnitudes. 
We made a cut in absolute magnitude because we did not 
have time to observe all 35 sources and by just selecting 
the most luminous we enable the luminosity distribution 
to be more similar to that of the control sample. Re- 
ichard et al. (2003) found that the optical colors of BAL 
quasars are redder than those of non-BAL quasars. This 
can naturally be explained by dust-reddening. A com- 
parison of composite quasar spectra formed by matched 
BAL and non-BAL SDSS sub-samples shows the mean 
amount of reddening is E{B — V) — 0.023. We determine 
the dereddened absolute S-band magnitudes of the BAL 
quasars assuming a rest-frame color oi B — i* = -1-0.35 
(Schneider et al. 2002) and incorporating an offset of 
—0.2 magnitudes to account for reddening. The absolute 
magnitude cut to restrict the sample to 30 quasars is at 
Mb < —26.6. Only three of these 30 BAL quasars are 
detected in the radio. These three quasars have 1.4 GHz 
flux-densities in the range 2-10 mJy, so the contribution 
due to synchrotron emission at 850 fim is expected to 
be much less than 1 mJy and hence can be neglected 
(Willott et al. 2002). 

2.2. Observations 

The entire sample of 30 BAL quasars was observed 
in photometry mode with the SCUBA bolometer ar- 
ray at the JCMT from February to June 2003. Ob- 
servations were made simultaneously at 850 /xm and 450 
fim. The majority of the observations were performed in 
good to average conditions (atmospheric zenith opacity 
at 1300 fim in the range 0.05—0.12) and hence the 450 fim 
data arc not nearly as effective at detecting thermal dust 
emission as the 850 /im data. Henceforth we will concen- 
trate on the 850 /xm data, but the 450 fim data from some 
observations are also given in Table 1. The source was 
chopped between two bolometers on the array to increase 
the effective time on source. 

The data were reduced using the SURF package, fol- 
lowing mostly the same methods as in Archibald et al. 
(2001). The main difference is that the sky opacity was 
determined from polynomial fits to cither the JCMT wa- 
ter vapor monitor or, on the occasions when this was not 
operating, the NRAO 350 /um tau meter. The data for 
the two bolometers were reduced separately and the fi- 
nal fluxes measured with each bolometer were combined 
weighting by the inverse of the variance. For some obser- 
vations the secondary bolometer data arc unstable and 
for these only data from the primary bolometer is used. 

The resulting flux-densities are given in Table 1. We 
have achieved our aim of reaching an rms sensitivity at 
850 /im lower than 3.33 mJy (i.e. 3cr = lOmJy) for all of 
the 30 broad absorption line quasars observed. In Table 1 
we highlight in bold the eight quasars which are detected 
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Table 1. Table of observed objects. 



Source 


z 


Mb 


BI / kms-i 


Type 


5850 / mjy 


5450 / ™Jy 


SDSS J104109.85+001051.8 


2.250 


-26.68 


1913 


Hi 


0.54 + 2.26 




SDSS J104130.17+000118.8 


2.068 


-27.66 


672 


Lo? 


2.04 ± 2.73 




SDSS J104233.86+010206.3 


2.123 


-27.20 


401 


Hi 


-0.56 + 2.80 




SDSS J104841.02+000042.8 


2.022 


-26.74 


1176 


Hi 


4.22 + 2.00 




SDSS J110041. 19+003631. 9 


2.017 


-27.30 


4687 


Hi 


0.03 + 2.13 




SDSS J110623.52-004326.0 


2.450 


-27.80 


4034 


Lo? 


-0.93 + 3.26 




SDSS J121633.90+010732.8 


2.018 


-27.04 


2200 


Hi 


3.45 + 2.62 


26.7+14.9 


SDSS J121803.28+001236.8 


2.010 


-26.90 


269 


Hi 


5.44 + 2.47 




SDSS J122228.39-011011.0 


2.284 


-26.70 


678 


Lo? 


-0.32 + 2.15 




SDSS J123056.58-005306.3 


2.162 


-26.91 


652 


Hi 


1.29 + 3.05 




SDSS ,1124720.27-011343.1 


2.283 


-26.81 


336 


Hi 


6.88+1.91 




SDSS ,11318.53.45+002211.4 


2.079 


-26.97 


1473 


Hi 


4.45 + 3.08 




SDSS ,1134145.13-00,3631.0 


2.205 


-27.17 


870 


FcLo 


0.17 + 2.66 




SDSS ,1134544.55+002810.8 


2.516 


-27.46 


1510 


Hi 


-3.14 + 2.54 


1.1 + 11.8 


SDSS J135317.80-000501.3 


2.320 


-26.73 


9821 


Lo 


4.17 + 2.73 




SDSS J135559.04-002413.6 


2.332 


-27.65 


1525 


Hi 


4.15+1.94 




SDSS J135721. 77+005501.1 


2.001 


-27.48 


235 


Hi? 


4.14 + 2.60 


39.5 + 21.8 


SDSS J140918.72+004824.3 


2.000 


-27.02 


1411 


Hi 


9.99 + 3.01 


23.1 + 15.7 


SDSS J142050.34-002553.1 


2.103 


-26.95 


3442 


Lo? 


3.47 + 2.33 


-26.9 + 25.8 


SDSS J143022.47-002045.2 


2.544 


-26.61 


1957 


Hi 


-1.62 + 2.78 




SDSS J144256.86-004501.0 


2.226 


-27.51 


815 


Hi 


4.33 + 2.35 


9.6 + 21.7 


SDSS J145045.42-004400.3 


2.078 


-27.27 


238 


Hi 


1.85 + 2..38 


-31.8 + 23.1 


SDSS J150033.52+003353.7 


2.451 


-27.86 


4257 


Hi 


8.13 + 2.56 




SDSS J150206.66-003606.9 


2.202 


-27.17 


406 


Hi 


0.29 + 2.51 


11.9 + 22.2 


SDSS J151636.79+002940.4 


2.240 


-28.44 


4035 


Lo 


-1.37 + 2.50 


3.0+ 15.0 


SDSS J152913.85-001013.8 


2.073 


-27.31 


209 


Hi 


-1.31 + 2.52 


30.5 + 22.4 


SDSS J170056.85+602639.8 


2.125 


-27.08 


1400 


Hi 


6.62 + 2.64 




SDSS J171652.35+590200.2 


2.369 


-27.16 


620 


Hi 


10.53 + 3.04 




SDSS J172012.40+545601.0 


2.099 


-27.33 


1249 


Hi 


3.75 + 2.06 


27.9 + 21.9 


SDSS J172341.09+555340.5 


2.113 


-27.11 


3497 


FeLo 


-3.74 + 2.28 


-10.4+17.6 



Note. — SCUBA photometry of SDSS broad absorption line quasars at 2 < ^ < 2.6 drawn from the sample of Reichard et al. (2003). The 
S-band absolute magnitudes are derived from the SDSS EDR z*-band absolute magnitudes by (i) a correction to our adopted cosmology; (ii) a 
color correction of — z* — +0.35 (Schneider et al. 2002); (iii) dereddening by —0.2 magnitudes to account for the average amount of intrinsic dust 
reddening compared to non-BAL quasars. Balnicity indices BI are for the C IV line and were computed by fitting a composite quasar spectrum to 
the continuum (column 12 of Tabic 1 of Reichard ct al. 2003). Type denotes whether the quasar is a Hi-ionization, Lo-ionization or Fc-Lo-ionization 
BAL (column 18 of Tabic 1 of Reichard ct al. 2003 with '?' from column 17 of the same tabic). 850 /^m detections at the > la level arc highlighted 
in a bold font. The three quasars detected in the radio are SDSS J135721. 77+005501.1, SDSS J150206. 66-003606. 9 and SDSS J151636. 79+002940.4. 
None of them are detected at > 2(t at 850 jum. 450 /jm flux-densities are only given for those observations performed in good conditions where the 
data for the 450 ^ra array is reliable. There are no > 1(j detections at 450 jum. 



at the > 2fT level. Of these, four are securely detected 
at the > 3c7 level. A quick glance at the tabic shows 
that the detected quasars span a range of z, Mb and 
balnicity index. A full discussion of the sample statistics 
and a comparison with the control sample is pursued in 
the following section. None of the quasars are detected at 
even the > 2(t level at the shorter wavelength of 450 /im. 
This is not surprising given the sensitivity level reached. 
These data do not provide useful constraints on the sub- 
millimeter spectral indices. 

We estimate the submillimeter /far-infrared luminosity 
from the sub-millimctcr fluxes by assuming the isother- 
mal dust template of Priddey & McMahon (2001). This 
template consists of isothermal optically-thin dust with 
temperature T = 41 K and emissivity index /3 = 1.95. 
This gives a conversion between Lfir and S%^q of Xfir ~ 
2 X 10^^ (S'850 / mJy) ^0 at rcdshifts z > 1. However, 
as discussed in Willott ct al. (2002) these values of 
.^FiR, should be treated as lower limits since isothermal 
models likely underestimate the total far-infrared lumi- 
nosity. The four BAL quasars we have detected with 
SCUBA at > 3(7 signiflcancc have far-infrared luminosi- 
ties > 10^"^ Lq. In the rest of this paper we will use the 
observed submillimeter flux-density in preference to the 
far-infrared luminosity and thereby avoid uncertainty in 



Z/FiR due to the unknown dust spectral shape. 

3. DO QUASARS WITH BALS DIFFER FROM THOSE 
WITHOUT? 

3.1. Control sample 

In this section we combine our new data with data from 
the literature to determine how the sub-millimeter emis- 
sion depends upon whether a quasar shows broad absorp- 
tion lines or not. As mentioned previously, our sample of 
BAL quasars was defined so as to be as similar as possi- 
ble in terms of redshift and luminosity as quasars already 
observed in the sub-millimeter or millimeter. Priddey et 
al. (2003) obtained 850 /im data for 57 quasars in the 
redshift range 1.5 < z < 3 with SCUBA on the JCMT. 
Their sample is drawn from the Large Bright Quasar Sur- 
vey (LBQS; Hewett, Foltz & Chaffee 1995) and Hamburg 
Quasar Survey (HS; Hagen et al. 1999). Omont et al. 
(2003) obtained 1.2 mm data for 35 quasars in the red- 
shift range 1.8 < z < 2.8 with MAMBO at the IRAM 
30-m telescope. Their sample was drawn from the; Veron- 
Cetty & Veron (2001) catalog and has some overlap with 
the HS quasars observed by Priddey et al. The 1.2 mm 
fluxes and errors from Omont et al. were multiplied by 
a factor of 2.5 to estimate the corresponding values at 
850 yum (Omont et al. 2003; see also Eales et al. 2003). 
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In cases where both papers present data on the same 
source, we take the data from Priddey ct al. in favor 
of that of Omont et al., since these observations were 
made at the same wavelength as our data and do not 
require multiplication by a flux correction factor. Note 
that there are a couple of cases of detections by Omont 
ct al. and non-detections by Priddey et al. However, 
to take the Omont et al. data for these sources would 
bias the sample unfairly. Both these papers contained 
a few sources showing BALs which we exclude, since we 
require a non-BAL control sample. Following Priddey et 
al. (2003) and matching to the sensitivity reached for 
the BAL sample, we exclude from further analysis the 
three sources from the Priddey et al. sample and the 
two sources from the Omont et al. sample for which the 
observations were noisy (rms at 850 /i^m > 3.33 mJy). In 
summary, our control sample at 2 < z < 2.6 contains 25 
quasars from Priddey et al. and 16 quasars from Omont 
et al. 

Unfortunately for us, both the above surveys concen- 
trated on the most optically-luminous quasars to enable 
comparison with high-luminosity quasars at z ^ A. Our 
sample is drawn from the SDSS which has fainter lim- 
iting magnitudes than the above surveys and hence our 
sample contains typically less-luminous quasars. There 
has been no systematic sub-millimeter /millimeter survey 
of z » 2 SDSS quasars which would have provided the 
ideal control sample to ours. The median B-band ab- 
solute magnitudes for the samples are Mb = —27.2 for 
the BAL sample and Mb = —28.2 for the non-BAL sam- 
ple. Applying the Kolmogorov-Smirnov statistic to the 
two samples, wc find that their distributions in Mb dif- 
fer very significantly (probability 10^^%). We can also 
use the Kolmogorov-Smirnov statistic to determine the 
offset in Mb to be applied to the BAL quasar sample 
to make its Mb not signifcantly different from the non- 
BAL sample at the 5% level. We find that a shift of 
5Mb — —0.8 is needed to make the two samples have 
similar flux distributions. 

Omont et al. (2003) combined several different samples 
to investigate whether there is a correlation between the 
optical and sub-millimctcr luminosities of quasars and 
found no statistically significant correlation between the 
two luminosities. However, they also noted that at z « 
4 the weighted mean sub-millimeter flux of an SDSS- 
based sample (Carilli et al. 2001) was half that of the 
weighted mean sub-millimcter flux of a more luminous 
sample (median Mb brighter by 1 magnitude; Omont 
et al. 2001). This is suggestive that there may be a 
weak correlation between the optical and sub-millimeter 
luminosities that is shrouded by the large scatter. We 
will return to this issue later. 

We have also checked to see if the two samples arc 
matched in terms of redshift. The rcdshift ranges for the 
two samples are identical with 2.0 < z < 2.6, however 
the BAL sample has a median z = 2.16 and the non-BAL 
sample a median z = 2.27. Applying the Kolmogorov- 
Smirnov statistic to the two samples, we find that the 
probability that the two samples arc drawn from the 
same distribution in redshift is 4%. There is evidence 
here that the samples are not ideally matched in red- 
shift, but the difference in cosmic time between the me- 
dian redshifts is minimal and there can not have been 
significant evolution of the universe over this period, so 




-10 r 

-26.5 -27.0 -27.5 -28.0 -28.5 -29.0 -29.5 -30.0 

Fig. 1. — 850 /^m flux-density against absolute magnitude for 
the BAL quasar sample (flUed circles) and the non-BAL control 
sample (open triangles). Error bars show the Icr uncertainty on 
the measured fluxes. A comparison between the two samples and 
possible correlations in this plot are discussed in Sec. 3. The dotted 
line shows a positive correlation between sub-millimeter flux and 
optical luminosity where a diflference in Mb of 1 equals a change 
in Ssso of a factor of 1.5. 



this will not affect our results. 

3.2. Statistical analysis 

The main problem we are interested in is whether 
there is a difference in the sub-millimeter properties of 
BAL and non-BAL quasars. We have used the observed 
850 /im fluxes to calculate the mean, standard error on 
the mean and median fluxes for four different samples: 
(i) the BAL quasars of this paper; (ii) the combined Prid- 
dey et al. and Omont et al. sample of non-BALs in the 
same redshift range; (iii) only the Priddey ct al. non- 
BALs; (iv) only the Omont ct al. non-BALs. We have 
also used the observed 850 ^m fluxes with their statis- 
tical uncertainties to calculate the weighted mean and 
its associated error. These results are presented in Ta- 
ble 2. The median sensitivity reached in our survey is 
similar to that in the non-BAL control sample (P-l-0). 
Both the mean fluxes and the weighted mean fluxes of 
the two samples are comparable within the uncertain- 
ties. The non-BAL sample actually shows slightly higher 
mean fluxes although the BAL sample has a higher me- 
dian flux. Also shown in the table are the percentages of 
sources which are detected at > lOmJy and > 6.66 mJy. 
We do not consider the percentage of > 3a and > 2a 
detections since these numbers depend critically on the 
sensitivity distributions of the samples. We find that the 
BAL sample contains a lower percentage of these sub- mm 
luminous sources than the non-BAL sample. Howc;ver, 
the percentage is similar to that seen in the Priddey et 
al. sub-sample. 

Due to the large number of 850 /xm upper limits in 
both samples, we also use survival analysis statistical 
tests which can account for these limits (Feigelson & 
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Table 2. Statistics of quasar samples at 2 < z < 2.6 surveyed in the 
(sub)-millimeter. 



Reference QSO type A*" median rms < Sggo > Weighted < 5g5o > median Sgso %{S850 > 10) %('S'g5o > 6.66) 

This paper BAL 30 2.54 2.56 ± 0.67 2.55 ± 0.45 3.5 3 13 

P+O Non-BAL 41 2.70 3.34 ± 0.64 3.61 ± 0.37 2.3 7 20 

P Non-BAL 25 2.80 2.88 ± 0.79 2.83 ± 0.55 2.2 4 16 

O Non-BAL 16 2.25 4.05 ± 1.09 4.29 ± 0.51 3.5 13 25 



Note. — Sub-millimeter properties of various samples of quasars at rcdshift 2 < 2 < 2.6. All fluxes iu this table, Ssso, are at 850 /^m and are 
in units of mjy. References: P - Priddey et al. (2003); O — Omont et al. (2003); the subsets of the data in these papers used here is described in 
Sec. 3.1. A'^ is the number of the sources in the sample used for the statistical analyses in this paper. < Sgso > is the mean observed sub-mm flux 
of the sample (not weighted by the noise) and the associated standard error. Weighted < Ssso > is the weighted mean sub-mm flux calculated as 
Yl{wiSi) /T:{'Wi) where the weights arc w — 1 / and the uncertainty in < 5's5o > is [J2{wi)]^'^'^ . 



Nelson 1985; Isobe, Feigelson & Nelson 1986). To min- 
imise the number of upper limits, we consider sources 
with snr> 2 to be detections for these tests. We wish 
to determine whether the 850 /xm flux distributions of 
the BAL and non-BAL quasars are different. Using a 
variety of tests (the Gehan, logrank and Pcto-Prentice 
tests), the returned probabilities that the fluxes of the 
BAL and non-BAL samples are drawn from the same 
distribution lie in the range (54-99%). Again, we find no 
evidence that BAL and non-BAL quasars differ in their 
sub-millimeter fluxes. 

3.3. Effects of a correlation between sub-millimeter and 
optical luminosities 

In the previous section we performed several different 
statistical tests and found with each one that there is no 
difference in the sub-millimeter properties of our sam- 
ples of BAL and non-BAL quasars. As mentioned in 
Sec. 3.1, the fact that the luminosity distributions of the 
two samples are not exactly matched leads to the un- 
certainty that a positive correlation between optical and 
sub-millimeter luminosity combined with enhanced sub- 
millimeter luminosity in BAL quasars could have led to 
the observed results. We now consider whether a correla- 
tion between the optical and sub-millimeter luminosities 
exists and, if it does, whether it could have affected our 
results. 

We use the survival analysis bivariate correlation and 
regression tests (Isobe, Feigelson & Nelson 1986) to as- 
sess the signficance of any correlations and their slope in 
these datasets. For the sample of 30 BAL quasars, the 
probability that there is no correlation between 850 fxni 
flux and Mb is 45% for Cox's proportional hazard model 
and 72% for the generalized Kendall's tau test. The best- 
fitting slope for the regression of logj^g S'sso and Mb is 
0.03 ± 0.12. Therefore, there is no correlation present 
between these quantities for the BAL sample. 

We now check for a correlation between these same 
variables in the non-BAL sample formed from the com- 
bination of the Priddey et al. (2003) and Omont et al. 
(2003) data. We now find a very significant (probabilities 
of 0.05% and 0.3% using the two tests) anti-correlation 
between these quantities, i.e. a positive correlation be- 
tween the sub-millimeter and optical luminosities. This 
is quite surprising since Priddey et al. and Omont et 
al. performed similar statistical tests on each of their 
full samples and did not find any significant correlation. 
The reason that we have now detected this correlation 



may be that by combining the samples over a restricted 
redshift range we have disentangled luminosity and red- 
shift effects. 

It is important to remember that some of the most 
luminous observed quasars at these redshifts are grav- 
itationally lensed and if a substantial number of these 
quasars are magnified by even modest factors of « 2, 
then this would induce a correlation between Sgso and 
AIb even if one does not exist intrinsically. We do not 
know which of these quasars are gravitationally lensed, 
but repeat the correlation tests for subsets of the data 
in an attempt to minimise the effects of lensing. First, 
we repeat the tests excluding the two outliers on Fig. 

the source with 6*850 — 17.2 mJy and the one with 
AIb — —29.6. We still find significant probabilities of 
0.9% and 2% that there is no correlation present. Lin- 
ear regression with the parametric EM algorithm gives 
the slope of the correlation between log^o Ss5o and Mb 
as —0.23 ± 0.10. This means that a difference in Mb of 
1 equals a change in Ss5a of a factor of 1.7 ± 0.4. We 
now repeat the tests using just the 30 non-BAL quasars 
with the faintest absolute magnitudes, since if lensing is 
an issue it is most likely to affect the quasars with the 
brightest observed absolute magnitudes. The correlation 
probabilities are now 10% and 14% for the two tests - 
suggestive of a correlation but not significant at the 5% 
level. The slope of the correlation between logj^Q 6*850 and 
Mb is —0.16 ±0.17, i.e. a difference in Mb of 1 equals a 
change in Ss5o of a factor of 1.5 ± 0.5. This correlation 
is shown as a dotted line in Fig. ^ 

We conclude that whilst it is highly likely that a 
correlation between 5*850 and Mb does exist at some 
level, the scatter and weakness of the correlation mean 
that it is hard to prove significantly with current data 
for optically-selected quasars. Similar conclusions were 
reached by Omont et al. and we refer the reader to 
their discussion on the likely causes of such a correlation. 
Grimes et al. (2003) have investigated the Ss5q-Mb cor- 
relation for radio- loud quasars, which assuming that 5*350 
is not strongly influenced by whether or not a quasar has 
radio jets, would be indicative of the relation for quasars 
in general. They find that 5350 drops by a factor of 1.5 
per unit drop in Mb- This kind of correlation is similar 
to that found in the non-BAL sample. 

At the low-luminosity end of the quasar luminosity 
function, sources are generally too weak to be detected 
individually with an instrument such as SCUBA. Stack- 
ing analyses of faint, soft X-ray sources from Chan- 
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dra and XMM-Newton show a positive mean flux. Com- 
bining the results of Barger et al. (2001) and Waskett et 
al. (2003) gives a mean sub-millimeter flux for faint soft 
X-ray sources of 5*850 = 0.60 ± 0.15 mJy. The optically- 
unobscured AGN in these samples have Mb ~ — 21. This 
is 7 magnitudes fainter than the non-BAL sample which 
has a mean sub-millimeter flux of 5*850 = 3.34 ± 0.64 mJy. 
The correlation between 585o and Mb implied by these 
two points is that 585o drops by a factor of 1.3 per unit 
drop in Mb- This is somewhat less strong than the fac- 
tor of 1.5 per unit drop in Mb found at the brighter end 
of the luminosity function, but is further evidence that 
the correlation is likely not steeper than this. 

For the purposes of our BAL quasar survey we wish to 
know if the existence of a correlation between 5350 and 
Mb could affect our findings. The median absolute mag- 
nitude of the BAL sample is Mb = —27.2 and for the 
non-BAL sample it is Mb = —28.2. Therefore, we scale 
the observed sub-millimeter fluxes of the BAL sample by 
a factor of 1.5 (except for those with negative flux) to 
simulate the effects of a correlation with slope similar to 
that observed in the non-BAL sample. We repeated the 
two sample tests with these scaled BAL fluxes and the 
comparison sample of non-BAL quasars. The returned 
probabilities that the sub-millimeter fluxes of the BAL 
and non-BAL samples are drawn from the same distri- 
bution are 14% and 27%. This test shows that, even 
after allowing for the effects of a weak correlation be- 
tween sub-mm flux and Mb, there is not a significant 
difference between BAL and non-BAL quasars. 

3.4. Correlation with balnicity 

The new sample of SDSS BAL quasars we are using 
was selected using methods that are very sensitive to the 
presence of weak BALs (Reichard et al. 2003). Quasars 
with extremely weak BAL features (BI < 200kms~^) 
were excluded from our sample, but 43% of our sam- 
ple is made up of quasars with moderate balnicities of 
200 < BI < 1000 km s"^. It is therefore important to 
check whether our results could be affected by a corre- 
lation between sub-millimeter emission and the balnic- 
ity index. In Fig. |21 we plot the 850 fim flux against 
the balnicity index for our sample. The control sam- 
ple is not shown since these are non-BAL quasars with 
BI = Okms^^ and lie to the left of the plotted re- 
gion. It is clear that there is no strong correlation be- 
tween 850 fim flux and BI and both the four > 3a de- 
tections and the sources with negative measured fluxes 
are spread throughout the range of BI. Using the sur- 
vival analysis correlation routines we find that both the 
Cox and Kendall's Tau tests show there is no correla- 
tion between BI and 850 fim flux (the two tests give 
54% and 72% probabilities for no correlation present, re- 
spectively). One can also divide the two samples at the 
median balnicity of BI = 1400 km s~^ and calculate the 
weighted means of the low- and high-BI sub-samples as 
2.99 ± 0.62 mJy and 2.07 ± 0.65 mJy, respectively. Once 
again, there is no indication of a correlation between BI 
and 850 fim flux. 

BAL quasars are usually classified into Hi- and Lo- 
ionization BALs, dependent upon whether low- ionization 
broad absorption lines such as Mgii and Aim are 
present. LoBALs are more prevalent among infrared- 
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Fig. 2. — 850 /im flux-density against balnicity index BI for the 
BAL quasar sample. There is no correlation present between BI 
and sub-millimeter flux (details are given in Sec. 3.4). 



selected samples (Boroson & Meyers 1992) and are sub- 
ject to greater dust reddening than HiBALs (Reichard et 
al. 2003), implying they could be enshrouded in a dusty 
cocoon during an early evolutionary stage of quasar ac- 
tivity (Voit, Weymann & Korista 1993). It is therefore 
interesting to see if the sub- millimeter properties are cor- 
related with the classification. In Table 1 we give the 
BAL classifications for these quasars from Reichard et al. 
There are 22 HiBALs, 6 LoBALs (of which 4 have uncer- 
tain low-ionization troughs and definite high-ionization 
troughs) and 2 FeLoBALs (see Hall et al. 2002). All of 
the eight > 2a detections are in the 22 HiBALs. This 
is the opposite result to that expected if LoBALs re- 
side in more dusty environments. However, given the 
uncertainty surrounding the classiflcation of some of the 
LoBALs and the fact that Mgii moves out of the SDSS 
spectral range at z > 2.26 means that these classifica- 
tions are not too reliable and no firm conclusions can be 
drawn. It is certainly possible that if a sample of LoB- 
ALs were observed in the future with SCUBA that they 
could have a different distribution of sub-millimeter flux 
than non-BAL quasars. 

4. INTERPRETATION 

The analyses presented in Sec. 3 have clearly shown 
that quasars with broad absorption lines show the same 
sub-millimeter properties as quasars without broad ab- 
sorption lines. What does this tell us about the nature of 
the BALs and the evolution of quasars? The first ques- 
tion we are attempting to answer is "are BAL quasars 
linked to non-BAL quasars by an evolutionary sequence 
or by orientation?" . Since sub-millimeter emission is 
optically-thin there should be no orientation dependence 
of sub-millimctcr emission. Our finding that the sub- 
millimeter emission is not related to the presence of a 
BAL is therefore consistent with the orientation hypoth- 
esis that all quasars contains BALs, but only a fraction 
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of them are visible along our line-of-sight. 

The evolutionary hypothesis is slightly more compli- 
cated since there are no reliable models for the relative 
timings of star-formation and AGN activity and the as- 
sociated sub-millimeter emission and BAL formation. In 
the evolutionary scenario, broad absorption lines are ob- 
served at a special epoch in quasar evolution when two 
conditions are met: (i) the quasar luminosity becomes 
high enough to produce powerful winds that are capa- 
ble of accelerating gas to several thousand kms~^; (ii) 
the nuclear regions contain a considerable amount of dif- 
fuse gas and dust. The black holes and hence maximum 
attainable luminosities of quasars grow with time as mat- 
ter is accreted. Only when the luminosity reaches a cer- 
tain value will the quasar be capable of producing strong 
winds. One can imagine a scenario where there is a brief 
period after this luminosity has been attained in which 
the bulk of the diffuse gas is ejected from the galaxy and 
the object shows broad absorption lines along most lines- 
of-sight (see Fabian 1999 for rough calculations). After 
this, the quasar will continue to shine for as long as it 
has fuel (which presumably is not ejected because it is 
in a denser accretion disc geometry) but no longer shows 
broad absorption lines. 

How does this picture relate to the expected sub- 
millimeter emission? Sub-millimeter emission comes 
from cool dust 50 K) which is heated cither by hot 
young stars or the active nucleus. The heating source for 
quasars with luminous sub-millimeter emission is a long- 
disputed issue (e.g. Rowan- Robinson 2000) but the mere 
presence of a large dust mass implies that star- formation 
did not cease a long time ago (Hughes, Dunlop & Rawl- 
ings 1997). Furthermore, the existence of large molecu- 
lar gas masses in some sub-millimeter-luminous quasars 
(e.g. Combes 2001) shows that these host galaxies are 
at early stages in their evolution. For our purposes it 
is enough to know that diffuse, widely distributed dust 
close to the nucleus and/or star- forming regions will pro- 
duce large amounts of sub-millimeter emission. The gas 
responsible for BALs also contains dust (the redder op- 
tical colours of BALs compared to non-BAL quasars is 
evidence of this; Reichard et al. 2003; Richards et al. 
2003). In the evolutionary scenario, the dust gets driven 
away from the galaxy along with the BAL gas and there- 
fore one would expect that BAL quasars should show 
considerably stronger sub-millimeter emission than non- 
BAL quasars. Our data do not support this simple model 
and show that the dust mass heated by star-formation 
and/or the AGN is similar for ByVLs and non-BALs. 

This does not mean that our results disprove the evo- 
lutionary scenario and support the unification by orien- 
tation of BAL and non-BAL quasars, merely that the 
particular picture outlined above is ruled out. One pos- 
sible modification is that the bulk of the dust which 
produces a sub-millimeter luminous phase is not co- 
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incident (in space and/or time) with the BAL gas which 
is ejected. For example, Archibald et al. (2002) present 
a model whereby the bulk of the star-formation and sub- 
millimeter emission from a forming galaxy occurs be- 
fore the black hole has grown large enough to become a 
quasar. Therefore when most quasars are shining in the 
optical they have already exhausted much of their dust 
which has been expelled by starburst winds or become 
locked up in stars. In this model, the BALs could still 
be formed at the special epoch outlined previously, but 
the dust mass at this epoch may not be that different to 
later on when a non-BAL quasar is observed. The main 
factor for whether a quasar shows strong sub-millimeter 
emission or not would just depend upon how de-coupled 
the star-formation and AGN epochs are. 

5. CONCLUSIONS 

We have conducted the first systematic survey of sub- 
millimeter emission from broad absorption line quasars. 

By comparison with matched samples of quasars which 
do not show the BAL phenomenon, we draw the following 
conclusions: 

• The sub-millimeter luminosities of BAL quasars 
are indistinguishable from those of quasars with- 
out BALs. 

• Since sub-millimeter emission is optically-thin and 
isotropic these results are consistent with the sce- 
nario where all quasars contain broad absorption 
line gas, but only a fraction with a particular ori- 
entation actually show BALs in their spectra. 

• The data are inconsistent with a model in which the 

BAL phenomenon indicates a special evolutionary 
stage which co-incides with a large dust mass in the 
host galaxy and a high sub-millimeter luminosity. 

• Our results provide constraints on alternative mod- 
els whereby BALs occur at a special stage in the 
evolution of quasars. 
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